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Despite decades of intensive studies, the effective interactions between strongly charged colloids 
still remain elusive. Here we show that a strongly charged surface with a layer of condensed counter- 
ions behaves effectively as a conductor, due to the mobile nature of the condensed ions. An external 
source charge in its vicinity is therefore attracted towards the surface, due to the image charge effect. 
This mechanism leads to correlational energies for counter-ions condensed on two distinct surfaces, 
as well as for free ions in the bulk. Generalizing Debye-Huckel theory and image charge methods, 
we analytically calculate these correlation energies for the two-plates problem, at the iso-electric 
point, where condensed counterions precisely balance the bare surface charges. At this point, the 
effective interaction between two plates is always attractive at small separation and repulsive at 
large separation. 

PACS numbers: 82.70.Dd, 83.80.Hj, 82.45.Gj, 52.25.Kn 
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The effective interaction between charged objects in- 
side an electrolyte or plasma is a problem of fundamen- 
tal importance to many branches of physical sciences. In 
past few decades, there has been accumulating experi- 
mental and numerical evidences [iHSl of effective attrac- 
tion between likely charged objects. This is in direct con- 
tradiction with the classical DLVO theory [10|, [ll[ , which 
always predict repulsion between two likely charged ob- 
jects. Among these experimental works, the recent work 
by Tata et. al. [9| is particularly striking: These authors 
directly imaged bound states formed by clusters of likely 
charged colloids in highly de-ionized water, with influ- 
ences of boundaries and multi- valence counter-ions care- 
fully excluded. The attraction appears to exist at scale 
of Debye length, which is several hundreds nanometers. 

Multiple mechanisms have been proposed to explain 
like charge attractions. Most of them arc based on 
the idea of counter-ion condensation, first discovered by 
Manning 1^ in highly charged lines. Shklovskii et. al. 



J13l Il4l | have shown that in the so-called strong cou- 
pling limit, where the Bjerrum length is much larger than 
the average spacing between condensed counter-ions, the 
latter form a strongly correlated 2D liquid, which lo- 
cally resembles a Wigncr crystal. Furthermore, the total 
charge of condensed ions may overtake the bare surface 
charges so that the effective charge density of a strongly 
charged surface may be reversed. This phenomenon is 
popularly called charge inversion or over-charging, and 
has been confirmed in numerical simulations. Further- 
more, previous studies ISMlSJ seem to indicate that in 
the strong coupling regime, correlations between densely 
packed counter-ions induce short range attraction be- 
tween two charged plates, with a characteristic scale of 
the Gouy-Chapman length, or inter-ion distance. 

Colloidal systems [3|, |j, |9| however are almost always 
in the opposite, weak coupling regime, and attractions 
were seen at much larger scale. Like-charge attraction in 



colloidal systems therefore can not be explained by the 
previous theories, and remains a fundamental challenge. 
The shear existence of like-charge attraction indicates 
that the colloidal systems must be close to the thresh- 
old of charge inversion. The mechanism of charge in- 
version in the weak coupling regime is presumably very 
different from that in the strong coupling limit, and may 
very well be of chemical nature. In the spirit of pragma- 
tism, we shall assume charge inversion docs happen in 
the weakly charged regime, and explore its consequences. 
The counter-ions close to the surface then form a dilute 
2d gas. We expect that a fraction of these ions can dif- 
fuse in the lateral directions. The free ions in the bulk, 
on the other hand, forms a dilute 3d plasma, and will be 
treated using standard Poisson-Boltzmann theory. The 
basic idea of this two-fluid model has been explored pre- 
viously by different groups |20l . |2l|, in various depths. 



By studying the responses of the 2d plasma to external 
perturbations, we can calculate the correlation energy of 
each ion due to its proximity to two distinct colloidal sur- 
faces, be it condensed on the surfaces, or freely diffusing 
in the bulk. This allows us to obtain analytic result about 
the correlation-induced effective interaction between two 
charged plates at the iso-electric point. 

The combination of a (negatively) charged colloidal 
surface and a layer of condensed counter-ions shall be 
called a dressed surface. The standard electrostatic 
boundary condition is 
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where $, $' are the potentials inside the electrolyte and 
inside the colloid, respectively, both evaluated near the 
surface. The dielectric constant of the colloid e' is typ- 
ically much smaller than that of water e w 80. It is 
therefore reasonable to take the limit e'/e — !• 0, so that 



$' drops out of the problem altogether. This approxima- 
tion however does not change our essential result. 

The net surface charge density i7net(2;) in Eq. ([1]) in- 
cludes both some fixed surface charges density and some 
mobile condensed counter-ions qnc{x). At the level 
of mean field theory, the latter is related to the lo- 
cal potential via the Boltzmann distribution: a^x) = 
gricCxp — /3q$(a;). Assuming weak perturbation, we may 
linearize this relation. At the so-called "iso-electric 
point" , the bare surface charges are completely neutral- 
ized by the condensed counter-ions, so that the dressed 
plate is overall charge neutral. The zero-th order term 
in (Tc(x) is then precisely cancelled by the fixed surface 
charge density. Therefore Eq. ([T]) is linearized into the 
following homogeneous Robin boundary condition: 
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If all condensed counter-ions were mobile in the trans- 
verse directions, the length scale /i would be just the 
Gouy-Chapman length, and takes the value of a few A 
to a few nm for strongly charged surfaces and therefore 
is much smaller than sizes of colloids and separation be- 
tween them. It is therefore attempting to ignore the left 
hand side of Eq. ([2]) altogether. The boundary condi- 
tion then reduces to that for a conductor. The idea that 
a charged surface with condensed counter-ions behaves 
as a conductor was first proposed by Shklovskii et. al. 
[ij, ll9[, in the context of charge inversion theory. It 
shall play important role in our analysis of counter-ion 
correlations below. 

To see how a dressed surface behaves as a conductor, 
consider a dielectric sphere (with radius a and e' ~ 0) at 
the iso-electric point inside a uniform external field i?o 
(no electrolyte). The potential outside the sphere can be 
expanded in terms of spherical harmonics: 
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where p is the induced dipole moment of the dressed 
sphere. Substituting this back into Eq. ^ we find: 
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For a typical colloid a ^ lOOnrn ^ fj, ^ A, hence the 
effective dielectric constant Coft of the dressed sphere is 
much larger than that of water. The sphere indeed be- 
haves like a conductor. 

To calculate the correlation energy for an ion, we shall 
follow Debye-Huckel's strategy. First consider an exter- 
nal point charge Q at a distance d from a dressed plate 
at the iso-electric point. As long as d 3> /i, we can again 
ignore LHS of Eq. ([2]) and treat the plate as a conduc- 
tor. Within linearized PB theory, then, the plate creates 



an opposite point image charge in the other side which 
exerts a reaction potential 
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on the source ion, and attracts it towards the plate. 

What if the source charge Q is one of the counter-ion 
condensed on the plate? In this case, the left hand side 
Eq. ([2]) can no longer be neglected. More importantly, 
the approximation of 2d fluid for condensed counter-ions 
is no longer valid at scale shorter than the average separa- 
tion between ions. Nevertheless, far away from the plate, 
the dressed plate should still behaves as a conductor. Its 
effect can be described by a point image —q, called pri- 
mary image, shown in Fig. [TJ The combination of source 
and images behaves effectively as a dipole, called "dipole 
zero" . In the regime of low density Uc^ a simple calcula- 
tion shows that the dipole moment p is given hy p ~ 2q^. 
If the density Uc is not low, however, the dipole moment 
p is rather difficult to calculate, and shall be reserved for 
a separate publication. At this stage, we shall treat p as 
a fitting parameter. 

Now introduce another identical dressed plates, i.e. the 
right plate in Fig. [l] at a distance d from the left plate. 
Within the linearized Poisson-Boltzmann theory, the ef- 
fect of the right plate is to introduce more image dipoles. 
To the lowest order, an image dipole at distance 2d to 
the right of the source arises, which shall be referred to 
as image dipole 1, see Fig. [TJ It has the same orientation 
and moment as the dipole 0. Image dipole 1 is reflected 
again by the left plate, and gives rise to the image dipole 
2, at a distance 2d to the left of the left plate. 

Dipole 1 and dipole 2 are responsible for the leading 
order correlation potential acting on the source charge 
q. Moreover, the interaction between the image dipole 2 
and the source charge q is identical to that between the 
image dipole 1 and the primary image. It follows that 
the lowest order correlation energy of the source charge 
is half [23[ of the interaction energy between dipole 
and the image dipole 1, which, in the linearized Poisson- 
Boltzmann theory, is given by 



2 Aire 



') 



=2rf 



87re4d3 



h:jTA^ + '^^d+2[ndf) 



-2Kd 



(6) 



As expected, this correlation energy is always negative. 
Since all counter-ions are identical, their total correlation 
energy is simply Eq. ^ multiplied by the total number 
of counter-ions on two plates. Furthermore it is rather 
straightforward to include all higher order image dipoles. 
The total correlation energy of condensed ions is: 



U, = - 4^ [Li3(e-2-'i) + 2KdLi2(e-2'''^) 
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FIG. 1: Two parallel plates, shown as vertical solid lines, are separated by distance d. A condensed ion q on the left plate is 
shown as a big black dot. (It is drawn away from the plate to make it more visible. ) It has a primary image —q (big open 
dot) just to the left of the left plate. The source charge and the primary image form the dipole 0, with a dipole moment p, 
pointing to the right. The right plate further creates an image dipole (dipole 1) with the same magnitude and orientation, at 
a distance 2d from dipole 0. The left plate further produces another image dipole (dipole 2), at a distance 2d to the left of the 
left plate. These two image dipoles are responsible for the lowest order reaction potential acting on the condensed ion q. The 
effects of higher order images dipoles are smaller by powers of e"^'"'. Also shown in the lower part are a free ion between two 
plates, as well as a few low order images. All images interact with source ion. 



where Li„j = J2kLi z'^/k"^ are polylog functions, and A 
is the total area of the plates. This result is in excel- 
lent agreement with direct Monte Carlo simulation of 
the implicit two fluids model, where the free ions be- 
tween two plates are treated using linearized Poisson- 
Boltzmann theory, as shown in Fig. [21 This clearly es- 
tablishes the image induced dipole-dipole interaction as 
a main mechanism of electrostatic correlation in strongly 
charged colloidal systems. 

The same method can be used to calculate the corre- 
lation energy of the free ions between two plates. A free 
ion (big closed circle in the lower part of Fig. [1]) between 
two plates has an infinite sequence of image charges in its 
each side. A few low order images are shown in Fig. [T] 
Since the dressed plates behave as conductors, all odd 
order images (small open circles) have opposite charges 
—q, and are at distances 221^2, 2(2;i_2 + d), 2(zi^2 + 2d), . . . 
from the source charge, where zi^2 are the distances from 
the source to the left/right plate respectively. All even 
order images (small closed circles) have like charges q, 
and are at distances 2d, 2 • 2d, 2 ■ 3d, . . . from the source 
charge. Summing up the interaction between the source 
charge and all its images, and integrating over the space 
between plates, we obtain the total correlation energy for 
all free ions between plates as [2J| 
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and therefore does not contribute to the effective inter- 
action between plates. The apparent logarithmic diver- 
gence at small z is due to our improper treatment of 
short scale details, but has no effect on our result. The 
second integral in Eq. (|8]) (due to even order images) can 
be easily calculated: 
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By simple variable transformations zi.2 + /d — > z, the 
first integral (due to odd order images) can be trans- 
formed into Jj, dz z~^ e~'^'^'^ , which is independent of d. 
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which, somewhat surprisingly, is always positive. Numer- 
ical simulation of this correlation energy shall be carried 
out in a forthcoming publication. 
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FIG. 2: Correlation energy of condensed counter-ions as a 
function of interplate distance. Dots: simulation of the im- 
plicit two- fluid model; line: analytic prediction Eq. ((7|. The 
data are renormalized to be unity at the shortest separation 
d — lOA. In this simulation, there are 40 point-like mono- 
valence counter-ions condensed on to each plate with size 
lOOA X iOOA. The boundary conditions are taken into ac- 
count using image charge methods. The system is deep in the 
weak coupling regime. 

The total correlation energy for the two-plates system 
is the sum of Eq. ([7]) and Eq. ^. At the level of our 



approximation, it is also the interacting free energy be- 
tween two charged plates. In the large separation regime 
nd > 1, t/c -- -e~'^'"^ / nd, while Uf -- e"^'"^. Therefore 
Uf always dominates and two plates repel each other. 
By contrast, in the small separation regime kc? <C 1, 
Uc ^ —l/nd^, while Uf ^ —\og{Kd). Hence Uc always 
dominates and two plates attract. The quantitative be- 
haviors of C/tot is controlled by one dimensionless param- 
eter a = ricP^ K^ / noq^ . Three representative cases of var- 
ious values of a arc shown in Fig. [3l 




FIG. 3: The total correlation energy for three cases. The unit 
of [/tot is q^noA/Ane. 

The result Eq. © describes the effective interaction 
between two neutral conductor plates. Interestingly 
enough, it also describes the interaction between two neu- 
tral dielectric plates with e' ^ e « 80. In the latter case, 
all images are likely charged, and the total correlation 
energy is still given by Eq. (|8]), but with the first integral 
changing sign. This integral is however independent of d 
as we have shown above. Therefore the total correlation 
energy is still given by Eq. (j9|). 

We emphasize that the essential ingredient of our anal- 
ysis is the fact that mobile condensed counter-ions make 
the charged surfaces resemble conductors. This is clearly 
independent of the dielectric constant of the charged col- 
loids themselves. Our approach therefore applies to ar- 
bitrary strongly charged objects in arbitrary solvent. 

In the regime of low density of mobile counter-ions, we 
have p ^ 2qfi oc l/ric, hence the prefacton of Eq. ([7]) 
is inversely proportional to mobile counter-ion density 
Uc- [25J Rather ironically, therefore, the correlation in- 
duced attraction is stronger when there are fewer mo- 
bile counter-ions. The surface charge density on a real 
surface should be heterogeneous at atomic scales. This 
may lead to pinning of most counter- ions at the threshold 
of iso-electric point, and therefore may greatly enhance 
correlation-induced attraction between two such charged 
colloids. 

We suspect that strongly charged colloids in highly 
deionized water may be very close to the iso-electric 
point, due to some unknown self-tuning mechanism, that 
is probably of chemical nature. It this is so, the correla- 
tional energies Eq. (O and Eq. (jH]) are still valid. Direct 



Coulomb interaction between two charged plates propor- 
tional to cr(^cr2^ should also be included, where 17(^2 ^^^ 
the renormalized surface charge density. Depending on 
the relative signs of ai^^2^ this could be either repulsive 
or attractive. For the case of colloids with finite sizes, 
additional attractions arise due to the overall polariza- 
tion of colloids by nearby charged colloids, whose mag- 
nitudes are proportional to ((7^2)^1 and is also doubly 
screened. This effect was analyzed by Levin [20| , as well 
as by Zhang and Shklovskii [22|. The quasi-conductor 
nature of the dressed surface also played an essential role 
in their analyses as well. We note that in the proximity 
of the iso-electric point, all these effects arc subdominant 
to the effective interaction derived in our work. 

Detailed understanding of short scale physics and 
chemistry is necessary in order to determine the values 
of parameters used in our two-fluids model, such as the 
cffective/renormalized charge density of a highly charged 
surface, as well as the effective density of mobile con- 
densed counter-ions. This demands close collaborations 
between theorists and experimentalists, as well as be- 
tween physicists and chemists. 
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The extra 1/2 is common to all the image charge inter- 
actions. 

Note that we assume that the correlation energy for each 
ion is sufficiently small that it does not change the uni- 
form density distribution of free ions between plates. 
This relation of course does not hold for arbitrary small 
Tic. Because if the length scale ^ in Eq. ((2]) becomes com- 
parable with the separation d, our dipole approximation 
of the image charges is no longer valid. 



